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Motivated by recent experiments [Zhang et ah, Phys. Rev. Lett. 95, 216801 (2005) and Ellenberger et at, 
cond-mat/0602271 1 reporting novel ringlike structures in the density-magnetic-field (n2D-B) diagrams of 
the longitudinal resistivity p xx of quantum wells with two subbands, we investigate theoretically here the 
magneto-transport properties of these quantum-Hall systems. We determine p xx via both the Hartree and the 
Kohn-Sham self-consistent schemes plus the Kubo formula. While the Hartree calculation yields diamond- 
shaped structures in the n^o-B diagram, the calculation including exchange and correlation effects (Kohn- 
Sham) more closely reproduces the ringlike structures in the experiments. 
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1 Introduction 

Quantum Hall systems, i.e., two-dimensional electron gases (2DEGs) under strong magnetic fields, exhibit 
a variety of novel physical phenomena [ 1 1. Besides the widely known quantum Hall effects, these systems 
also display quantum Hall ferromagnetism - an spontaneously formed spin-polarized many-body state 
arising from the interplay of the Zeeman, Coulomb, and thermal energy scales within the highly-degenerate 
Landau states of the 2DEG 0. 

Quantum Hall ferromagnetism has been experimentally investigated via magneto-transport measure- 
ments of the longitudinal p xx and Hall p xy resistances in both the integer and fractional quantum Hall 
regimes. A crucial ingredient in these experiments is the occurrence of opposite-spin Landau level cross- 
ings [3] near the Fermi energy. At these crossings the system may find it energetically favorable to spin 
polarize thus reducing its energy due to the Pauli-Coulomb exchange interaction. Resistance spikes (or 
"anomalous" peaks) are found in p xx at these crossings. Corresponding dips or bumps in p xy are some- 
times reported. Most importantly, these features in p xx and p xy are hysteretic in many cases thus signaling 
Ising-like ferromagnetism [4, 5, 6 , [3. 

Recent experiments in two-subband quantum Hall systems [8, 9| have found interesting ringlike struc- 
tures (Fig. 1(a), reproduced from Ref. [8|) in the density-magnetic-field (i^d-B) diagram of p xx . The 
authors of Ref. [ 8 1 claim that these unusual structures arise from quantum Hall ferromagnet states in the 
system, while those of Ref. [9] explain their ringlike structures in terms of a phenomenological single- 
particle picture. In both experiments crossings of Landau levels from distinct confined subbands occur at 
the Fermi level. 

Here we theoretically investigate the magneto-transport properties of quantum-Hall systems with two 
occupied subbands. We consider a modulation-doped GaAs/AlGaAs quantum well, similar to that of Ref. 
ISl. flOl . We determine the subband structure of the well via both (i) the Hartree approach and (ii) the 
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Experiment 
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Fig. 1 Experimental (a) (adapted from Ref. |8|) and calculated density-magnetic-field p xx diagrams within the 
Hartree approximation (b) and the SDFT/LSDA (c). Only the ringlike structure corresponding to filling factor v = 6 
is shown in (b) and (c). 



Spin Density Functional Theory (SDFT) ifTTI [within the Local Spin density Approximation (LSDA)] 
implemented via the usual Kohn-Sham scheme 1121 . We obtain the magneto-conductivity tensor (and its 
inverse, the magneto-resistivity) from the Kubo formula. We find diamond-shaped structures [Fig. 1(b)] in 
the U2D-B diagram within the Hartree approximation, while the SDFT-LSDA calculation yields smoother 
structures [Fig. 1(c)] more closely resembling the ringlike shapes in Refs. |8j[9l [e.g., Fig. 1(a)]. We do not 
find any hysteresis or discontinuities in the DFT-LSDA Landau-level fan diagram, which would indicate 
easy-axis (Ising) quantum Hall ferromagnetism in the system fl4j[5][6][7). However, we can not rule out the 
possibility for easy-plane ferromagnetism. 



2 Self-consistent approach 

Hartree and SDFT/LSDA calculations. Within the usual implementation of the SDFT/LSDA formulation 
in the context of the effective mass approximation, we self-consistently solve the Kohn-Sham equations 
for the conduction electrons in an external potential. Because of the translational symmetry in the xy plane 
of the 2DEG, the problem is separable [7 ]. For the z direction we have two (spin up and down) coupled 
one-dimensional Schrodinger equations 



h 2 d 

2m dz 2 



+ v e}f 0; i n h n l\) 



X? (z) = sfxf (z) , (1) 



with m being the effective mass, i = 1,2... the subband index, and v^Az) the effective single-particle 
potential. While the heterostructure potential confines the longitudinal motion, the perpendicular magnetic 
field collapses the transversal motion into Landau levels with energies e n = (n + 1/2)Tllu c , uj c = eB/m, 
each of which with a macroscopic degeneracy eB/h (~ 10 11 cm -2 ). Therefore, the Kohn-Sham single- 
particle energy of each electronic level is 

+ M huj c + y9VbB . (2) 

where gnBBa z /2 is the ordinary Zeeman term and g the effective Lande factor. The spin-dependent effec- 
tive potential in Eq. (fl~|i consists of two main contributions: the confining potential of the heterostructure 
v c (z) and the Coulomb potential which in the Kohn-Sham/SDFT scheme is split into Hartree Vh(z, [n]) 
and an exchange-correlation (XC) v xc (z, [nj, nj) contributions. The one-dimensional effective potential 
in the z-direction then reads 

v eff( z > l n h n l}) = v c(z) + v H (z, [n])+v£(z, [n T ,nJ) . (3) 
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In our simulations, the structural potential v c (z) represents the square quantum well investigated in Ref. J8) . 
The Hartree potential is calculated by solving the Poisson's equation taking into account the electronic 
charge density and the positive charge profile of the depleted modulation-doped regions. For the XC 
potential we use the Vosko-Wilk-Nusair (VWN) [ 1 3 1 parameterization of the LSDA. 

Since vjj and v xc are both functionals of n-p and n^, which are calculated from the Kohn-Sham orbitals 
(z), Eqs. (Q} and (f3]i are solved self-consistently. The densities, and any other observable of interest, 
are calculated using a Gaussian broadened density of states (DOS) weighed by the Fermi function (aufbau 
principle). Each level is characterized by a width T a ' with an electron density eB / h. Note that the Hartree 
calculation is obtained from the above self-consistent scheme by simply turning off the XC contribution. 

Magneto-transport. We determine the magneto-resistances from the magneto-conductivity tensor (p = 
er -1 ), calculated within the self-consistent Born-approximation model of Ando and Uemura lfl4l[T5l . For 
short-range scatterers, 
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where a x 



en 2 D/B 



Aa xy with Aa xy being a small correction |[T4l [151 . We model the extended 



Landau states by a Gaussian DOS <7 ex t(£) of width r oxt = 0.25 meV [ 16 1. Next we present our results. 



3 Results and discussions 

Parameters. We consider a GaAsZAlo.3Gao.7As modulation-doped quantum well HI with a width L = 
240 A and two symmetric <5-doped layers with planar density n^D = 10 12 cm -2 (which is identical to 
the electron density in our simulations), set back from the well by 240 A. We use an effective g-factor 
q = -0.44, m = 0.063, a potential well offset - 236.43 meV 031 [El and T = 340 mK in our 
simulations. 

Hartree calculation. First of all we note that in a Hartree-type calculation the Zeeman interaction is 
the only source of the spin splitting of the energy levels. Moreover, this splitting increases linearly with 
B and does not depend on the electron density n^D- Using standard parameters for GaAs wells, e.g. 
g = —0.44 and m = 0.063, our Hartree calculation yields 112D-B p xx diagrams showing small diamond- 
shaped structures (this shape follows from the linear dependence of the spin splitting with B) while the 
experiment [8 | exhibits much larger ring-like structures [Fig.QJa)]. In order to generate structures with an 
area comparable to the experimental rings, we need to use an enhanced g ~ —2. 1 in our Hartree simulation. 
This value is consistent with the experimental estimate of Ellenberger et al. [9| who correctly suggests (as 
our SDFT/LSDA calculation below shows) that this enhancement is due to exchange-correlation effects. 
Figure |TJb) shows the v = 6 diamond. The small area of the diamonds calculated with the standard 
parameters suggests that the Zeeman effect is not the only source of spin splitting in the experiments. This 
is indeed the case, as we find in the SDFT/LSDA calculation. 

SDFT/LSDA calculation. Here, in contrast to the Hartree case, the spin-dependent exchange-correlation 
potential (z, [rif, njj) in the Kohn-Sham equations gives rise to a density-dependent spin splitting, not 
linear in B. Figure[Tfc) shows the v = 6 structure in the n,2D-B diagram, calculated via the SDFT/LSDA 
with standard parameters for GaAs wells. We can see a smoother structure which more closely resembles 
the ring-like structures in the experiment of Ref. [8| [Fig. |TJa)]. We have also calculated the full range 
n,2D-B p xx diagram and find that it is very similar to the experimental one, Fig. [TJa) [ 19]. Therefore, the 
SDFT-type calculation with the standard GaAs parameters results in more rounded structures with areas 
and shapes similar to those reported in Refs. [8 1 and 

Quantum Hall ferromagnetism. The above results seem to corroborate the conjecture in Ref. [8 1 that the 
ringlike shapes of the structures should be due to density-dependent exchange-correlation effects. How- 
ever, our SDFT/LSDA simulations do not show ferromagnetic phase transitions. The Landau fan diagram 
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and the resistivities (p xx and p xy ) do not show any hysteresis or discontinuities which would signal an 
easy-axis (Ising-like) ferromagnetic transition ||4] [5] |6] |7) . In addition, the spin-polarization at the center 
of the rings are too small to be understood as an easy-axis ferromagnetic phase. However, we cannot rule 
out easy-plane (xy) ferromagnetism. In fact, an early experimental investigation |20| has reported on the 
observation of both easy-axis and easy-plane quantum-Hall ferromagnetism in similar systems. 

4 Summary 

We have calculated the U2D-B diagram of p xx within both the Hartree and the SDFT/LSDA approaches 
(plus the Kubo formula) in two-subband quantum Hall systems. While our Hartree calculation with an 
artificially enhanced g factor of —2.1 yields diamond-shaped structures, the calculation with density- and 
spin-dependent exchange-correlation effects and the standard g = —0.44 value, more closely reproduces 
the ringlike shapes seen experimentally. 
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